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Abstract 

The amplitude and transverse quadrature component squeezing of coherent light in high 
Q cavity by injection of atoms of two-photon transition are studied. 

The Golubev-Sokolov master equation and generating function approach are utilized to 
derive the exact variances of photon number and of transverse quadrature component as 
function of t. The correlation functions and power spectrums of photon number noise and 
of output photon current noise are also investigated. 


1 Introduction 

In this work, the amplitude squeezing as well as the transverse quadrature component squeezing of 
coherent light in high Q cavity by injection of atoms are investigated. The interaction is assumed 
to be two-photon transition type and the initial mean photon number N is assumed large. 

The interaction interval r for individual atom is taken the faverable value i r/ g, where g denotes 
the effective coupling constant between the atom and the single mode light. This value makes 
each incoming atom to emit two photons during passing the cavity. 

Our approach is based on Golubev-Sokolov master equation! 1 !. Since this equation was doubted 
by Benkert and Rzazewskb 2 ' for it may give negative probabilities. We will do some discussion on 
it first. To our view, even if Golubev-Sokolov equation does not have the meaning as a common 
differential equation, it is able to give correct mean values, variances and correlation functions of 
appropriate quantities, when it is utilized along with generating function method. In this work, 
this approach is not only used to derive photon number variance, the power spectrums of steady 
photon number noise and of output photon current noise, but also is generalized to study the 

squeezing of transverse quadrature component. 

In the investigation of photon number variance, we find that, in the case that the steady mean 
value of photon n 3 is much larger than the initial mean value of photon N, the ratio (A n(t) 2 )/ ( n(t )) 
will first drop to a value which is much smaller than its steady state value 1/2, and then turns up 
to approach 1/2. 
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In the investigation of squeezing of transverse quadrature component, we get that its variance 
square is expressed by ( n(t))/AN , hence the correspondent steady value njAN maybe either 
smaller on larger thatn the standard value 1/4. It is interesting to note, this steady value is 
related to the initial parameter N . Furthermore it does not depend on whether the injection is 
regular or poissonian. 


2 Model, Golubev-Sokolov master equation and gener- 
ating function approach 

We assume that the initial state of light in the high-Q cavity is single mode coherent light with 
mean photon number TV » 1. The injected two level atoms are in upper level, they interact with 
the cavity field by resonant two-photon transition: u 0 = 2u. 

The change of density matrix of photon field due to its interaction with a single atom initially 
in upper level is described by 


{up) rnn — p mn cos {gr\J (m + l)(m + 2)) cos (gr-\J (n + l)(n + 2)) 

+Prn- 2,n — 2 sin(pr \Jm{rn - 1)) sin(gT^Jn(n - 1)) - p mn . ( 1 ) 

For lage TV, m and n for important p mn are also large, so that gr^m + l)(m + 2) may be 
approximated' 3 ) by pr(m + ±), etc. If we take the value of r as J then eq.(l) turns out to be 

(up) mn = (— l) m "Pm— 2,n— 2 — Pmn (2) 

which means each atom emits two photons during passing the cavity, namely the quantum effi- 
ciency of photon production equals one. 

We assume that the atoms enter the cavity one by one and at most one atom in the cavity 
every moment. Therefore after injection of k atoms, p will change to (1 + u) k p. 

If the injection is of poissonian statistics with r as mean injection rate, the average number of 
injected atom during the interval t -*• t + dt will equal K = rAt. Thus' 11 


p(t + At) 


£e" K fp (1 +i)V = e rta V(i), 


(3) 


leading to 


dp(t) 

—JT = rup. ( 4 ) 

ill pump v / 

This is just the pumping term in the well known Scully-Lamb master quation. 

For regular pumping, k itself is a definite number, k = rAt with r denoting the injection 
rate.Therefore 


which leads to' 1 ) 


p(t + At) = (1 + u) rAt p(t) = e rAtln ^p(t), 


dp(t) 

dt 


I pump 


= r[fn(l + u)\p(t) 


(5) 

( 6 ) 
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for regular injection. By adding the cavity damping term, Golubev and Sokolov got the equation 
= r[ln(l + u)p(t)] mn + T[~(m + n)p mn {t) + ^(m + l)(n + l)p m+1>n+ i(i)], (7) 

in which F denote the cavity damping and the thermal photon is assumed negligable. 

Benkert and Rzazewski found 121 that this equation gives negative p nn when it is solved by 
letting f t p nn = 0 to derive the steady values of p nn . Let us see where this problem might come 
from. For regular injection rAf equals k therefore must be larger than 1. Thus A t cannot be 
taken as arbitrarily small. This in turn means eq.(7) may not be a differential equation of common 
sense, one ought to avoid by setting ^ be zero to get the steady value of p. Because of stepwise 
increase of rt , the strickly steady value of p may not exist. 

In practice, one usually only needs to calculate the expectation values, varicances or correlation 
functions of some relavant quantities. In this case it is better to evaluate these values directly 
rather than through evaluating p mn first. Generating function aproach is especially good for this 
purpose. In this work this approach will be used not noly to study the amplitude squeezing 
(photon number squeezing) but also generalized to study the squeezing of transverse quadrature 
component. 


3 Photon number squeezing^ 

Golubev and Sokolov, as well as some other authors, expanded the logurithm ln( 1 + u) and 
trancated at the second order of u: 

log(l +u) = u-^u 2 . (8) 

W.-h.Tan 151 and the present author 1 ® 1 has shown independently that for evaluating the variance 
square (A n 2 (£)), this treatment is correct. The result so obtained is identical to the exact solution, 
but it is not so for evaluating (A n(t) 3 ). In general, for calculating of (An(t) 1 }, one needs to expand 
ln( 1 + u) to l terms to get the correct value 1 ® 1 . 

As did in Ref[l], we introduce the generating function for (A n 2 (t)) as 

G(z,t) = ]Tp nn (t)z n } z <1. (9) 

n=0 

By utilization of eqs.(7) and (8), we get the equation for G(z,t) as 

= -1(3 - z 2 )(l - z*)G(z,t) + r(l - zf-^L. (10) 

This is a partial difference equation of first order, its general solution is expressed by one of 
its special solution mutiplied by the general solution of equation = T(1 — z ) 8< ^’^ . The 

latter will be determined by the requirement of initial condition. The desired solution so obtained 
expressed by 

G(z,t)=C(y, (11.1) 


391 



where 


y=l + (z- l)e- r ‘, 

/(*,<) = £[3(* - V) + \(z 3 - y 3 ) - i(z 3 - y 3 ) - i( 2 < - j, 4 )]. 
The values of ( n(t )) and (An 2 (f)) are easily obtained from G(z,t ): 


wt» = 


(An 2 (f)) = 


dz 
& G(z ,t) 
<9z 2 1 


.=1 + 


dGCM), 

dz 


U-i 


’ dG(z, t) 
dz 


1 2 


r /__ 2r\ _ r . r 

f +(jV- T )e + p. 


Z= 1 


-2rt 


The steady values of (A n(f)) and (An 2 (t)) exist. By letting t = oo, one gets 

(«)» = (An 2 (f)), = p. 


( 11 . 2 ) 

(11.3) 

( 12 . 1 ) 

( 12 . 2 ) 

(13) 


If we define 77(f) as (A n 2 (t)) / (n(t)) , then its steady value r] s will be the same as one photon- 
transition subpoissonian lasers. 

Eqs.(12) can be checked in the special case of ideal cavity (T = 0)! 3 1. 

The 7 ](t) defined above has different behavior for x(= nj/N) > 1 or < 1. In the latter case 
77(f) drops from its initial value and monotonically tends to the steady value 1/2. In the former 
case 77(f) first drops down to a minimum value T) min less than 1/2 and then turns up to approach 

1/2. For x » 1, T] m i n ~ yj C 1, therefore the correspondent state may be closed to the photon 
number eigen state. 

The steady state correlation function g(t) defined as following 


g(t) = tr [p s d// (O)a^ (f)a# (f)djy (0)] , f > 0 


(14) 


can also be evaluated by a generating function F(z,t). F(z,t ) satisfies the same differential 
equation as eq.(10), but has different initial conditions: 


F(7,0) = E("+l)Pn«.n+l- 2 " (15) 

n 

The g(t) so attained is 

0(0 = l«-i = ( n ) 2 s ~ * > 0. (16) 

The power spectrum of the steady state output photon current noise is related to g(t), in the 
case that the damping of the cavity field is mainly due to output, its expression will be 

PiH = ( 17 ) 

The correlation function (A7r(fi)An(f 2 )) for arbitrary t\ and f 2 can also be calculated by similar 
approach! 4 !. ^From it we obtain the power spectrum of steady state photonnumber noice as 

Pn{w) = ( n )s^j2 _|_ p 2 > (18) 

which mainly lies in low frequency region, in contrast to Pj(u ) given above. 
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4 The squeezing of transverse quadrature component. 

We are now generalizing the generating function approach to investigate the squeezing of quadra- 
ture components of a. 

Let a, the eigen value of a for the initial photon state, be real number, then a i = |(a + a*), 
0,2 = ^(a — a*) will be the longitudinal and tranverse quadrature components respectively. 

The mean value of longitudinal quadrature component is given by 

(ai(f)) = ^ ] \A; I 1 (0 • 


In our model (fip)n,n+i = — p n - 2 ,n-i — Pn,n-u which absolute value is not small as compared with 
|p nn+1 |. Actually it is almost twice as large as |p n , n+ i|. And the sign of [(1 + u) fc p]n,n+i varies 
alternately between positive and negative as k varies. Because of these features, the evolution of 
{ a.i(t )) could not be described by differential equations. The situation of transverse quadrature 
component is different. In our case { a 2 {t )) remains to be zero. We may generalize the generating 
function method to investigate its variance square, which is expressed by 

(A a 2 {t) 2 ) = \ + ^X^PnnW - ^S\/( n + 1 )(« + 2)p„ ) „ +2 (i)- (19) 

^ " n ^ n 


As before, £„ yj(n + l)(n + 2 )p nin+2 (t) may be approximated by £ n (n + §)Pn,n+ 2 (f)- Define 


G 2 (z,t) — y; Pn.n +2 (O^ 71 ) Z < 1, (20) 

n 

then 

x> + |w«(() = 2aAitl|„, + ?g 2 (z,oi„,. (21) 

We see that only first order derivative appears in eq.(21), therefore it is enoughl 6 ^ 4 ! to take just 
one term in the expansion of ln{ 1 +ff). The equation of G 2 (z, t ) can be derived accordingly, solving 
it as before, we get £ n (n + |)Pn.n+ 2 (^), which in turn yields 


(Aa 2 (t) 2 ) 


i n )s . _! /. _ ( n s) s _rt _ i n (t)) 
4JV 4 l N ’ AN ' 


( 22 ) 


This result may also be checked in the special case of ideal cavity. Setting e rt = 1 — Tt in eq.(22), 
we get 

= i + & 

which is the same as that given in Ref[3] by a completely different approach. 

One may show from eq.(22) that (A a 2 (t) 2 ) — \ may be positive or negative, depending on 
whether ( n s )/N is larger or smaller than one. The steady value of (A a 2 (t) 2 ) is given by 


(Aa 2 (t) 2 ) 4 


(n)s 
AN ’ 


(23) 
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which may be much larger or much smaller than 1/4. The latter means, in certain sense, deep 
suppression of phase noise. 

It is interesting to note that the stationary value {Aa, 2 (t) 2 ) s is still related to the initial pa- 
rameter N. 

It is also interesting to note that (Aa 2 (Z) 2 ), unlike (A n(t) 2 ), has no concern with whether the 
injection is regular or poissonian, since in the above derivation, Zn(l +u) is allowed to be replaced 
by u. 

This work is part of the project supported by the Chinese Doctoral Program Foundation of 
the Institution of Higher Education. 
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